A Hamiltonian approach for explosive percolation 
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We introduce a cluster growth process that provides a clear connection between equilibrium 
statistical mechanics and an explosive percolation model similar to the one recently proposed by 
Achlioptas et al. [Science 323, 1453 (2009)]. We show that the following two ingredients are essential 
for obtaining an abrupt (first-order) transition in the fraction of the system occupied by the largest 
cluster: (i) the size of all growing clusters should be kept approximately the same, and (ii) the 
inclusion of merging bonds (i.e., bonds connecting vertices in different clusters) should dominate 
with respect to the redundant bonds (i.e., bonds connecting vertices in the same cluster). Moreover, 
in the extreme limit where only merging bonds are present, a complete enumeration scheme based 
on tree-like graphs can be used to obtain an exact solution of our model that displays a first- 
order transition. Finally, the proposed mechanism can be viewed as a generalization of standard 
percolation that discloses an entirely new family of models with potential application in growth and 
fragmentation processes of real network systems. 

PACS numbers: 



The second-order critical point of percolation [1, 2] has 
been successfully used to describe a large variety of phe- 
nomena in Nature, including the sol-gel transition 
or incipient flow through porous media as well as 
epidemic spreading Q and network failure [1, 0, H, @- 
A long standing question of practical interest has been 
since, how the transition could be made more abrupt and 
in the limit become even of first-order. In other words, 
what ingredient must be tuned in the basic model of ran- 
dom percolation to change the order of the transition? 

Recently Achlioptas et al. 10] proposed a new mecha- 
nism on random graphs which they termed "explosive 
percolation" that exhibits first-order phase transition. 
Their process takes place in successive steps, with bonds 
being added to the system in accordance to a selection 
rule. At each step, a set of two unoccupied bonds are 
chosen randomly. From these two, only the one with min- 
imum weight becomes occupied. In Ref. [Toj ] . the weight 
is defined as the product of the sizes of the clusters con- 
nected by this bond (this is called "product rule"). Im- 
portantly, if the bond connects two sites that already 
belong to the same cluster, the weight is proportional to 
the square of the cluster size. Since unoccupied bonds 
connecting vertices in the largest cluster have the largest 
possible weight, these bonds will become occupied only 
if two of them are randomly chosen. Thus, this selec- 
tion rule hinders the inclusion of bonds connecting ver- 
tices that already belong to the largest cluster. As a 
consequence, bonds merging two smaller clusters will be 
selected more frequently, resulting in the fast growth ob- 
served. Their model was then implemented on a fully 
connected graph, however, it was shown that the same 
effect takes place on 2D square lattices [ll| as well as 
scale- free networks [HI, [l3[ • 

In this letter we investigate what are the basic princi- 
ples that lead to the first-order phase transition observed 
in the explosive percolation model. First we name merg- 
ing bonds those edges that connect vertices in distinct 




FIG. 1: Two ingredients for explosive percolation. Here we 
show a possible configuration for a growth process where, at 
each step, any unoccupied bond can be introduced in the 
graph. For instance, in this figure we show three bonds that 
could be added in the next step, namely, a, (3, and 7. The two 
ingredients for obtaining a sharp transition are the following: 
(i) bonds that keep the clusters approximately at the same 
size are favored over bonds that result in larger size discrep- 
ancies; and (ii) bonds that connect vertices in the distinct 
clusters (merging bonds) are favored over bonds that connect 
vertices in the same cluster (redundant bonds). Thus, among 
the bonds indicated, a has the smallest probability due to 
condition (ii), j3 is not accepted due to condition (i), and the 
most probable is the 7 bond. 



clusters, while redundant bonds are edges connecting ver- 
tices in the same cluster. We show that two conditions 
are necessary for obtaining a first-order transition in a 
growth process where bonds are included one by one, 
namely, the process has to favor the inclusion of bonds 
that keep all the clusters at about the same size, and the 
process has to preclude the introduction of redundant 
bonds, at least below the critical point. More precisely, 
merging bonds must be introduced with much higher 
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probability than redundant bonds. In Fig. [T] we show 
a pictorial description of these two ingredients. 

In order to validate our hypothesis, we propose an ex- 
tension of the percolation model that describes a gen- 
eral growth process in the space of graphs. For this, 
we define a Hamiltonian that depends on the graph G 
describing the network. The probability of finding the 
system in a certain state G will be given by P(G) = 
Z- 1 exp(-/3H(G)), where Z = £ G exp(-/3£T(G)). A 
simple form for a Hamiltonian that includes the two in- 
gredients is 

= (1) 

where the sum is over the entire set of clusters C, s, is 
the number of vertices in cluster i, and ti is the number 
of redundant bonds added to this cluster. If the number 
of bonds in the cluster is hi , we have = 1 + b$ — s, . Note 
that each time one includes a redundant bond, one also 
closes a new loop in the cluster, thus £j is also a measure 
for the number of loops in the cluster. 

We can now simulate a process of cluster growth con- 
trolled by the Hamiltonian of Eq. (fT]). This is performed 
by starting with a network of N vertices without bonds, 
so each vertex initially belongs to a different cluster. At 
each step, a new bond can be placed between any pair of 
vertices not yet connected. The probability of including 
a particular bond b is given by ~ exp(— pAHb), where 
AiJb is the energy change after including this bond. Such 
a growth model emulates equilibrium configurations of 
graphs following the Eq. |T]) and having a given num- 
ber of bonds Nb. However, since the removal/rewiring of 
bonds is not considered during growth, this corresponds 
to an out-of-equilibrium process. Consequently, some dif- 
ferences should be expected between the observed results 
and the actual thermal equilibrium. 

For small values of a, redundant bonds are favored over 
merging bonds, while for large values of a the opposite 
takes place. Let us investigate the asymptotic behavior 
in the two different scenarios. If redundant bonds are 
favored, one might expect that a new merging bond will 
be included only after the addition of all possible redun- 
dant bonds. Since clusters of equal size minimize Eq. (fT]), 
we can assume that, for low temperatures, all clusters 
have about the same size S, so that fully connected 
sub-graphs with S(S — l)/2 bonds are formed with I = 
(S— 2)/2. After adding the next bond, two of these 
clusters shall merge to form a new largest cluster, into 
which redundant bonds can be included. At this point 
we can calculate the energy variation for a redundant 
bond, AH r = (2S) a , and for a merging bond between 
pair of clusters, AH m = 2S 2 + (S - 1)(S - 2)(2 a - l)S a . 
Surprisingly, for any value of a, in the asymptotic limit of 
very large clusters, S — > oo, merging bonds have higher 
energy variation than redundant bonds, and the growth 
process with fully connected clusters is stable. 

The situation becomes quite different when the pres- 
ence of merging bonds is favored. As before, we use that 
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FIG. 2: Transition to explosive percolation. When the pro- 
cess favors redundant bonds, a — 2.5, the largest cluster fol- 
lows a slow continuous growth for the largest cluster. When 
the merging bonds are favored, a = 2.7, the system displays 
an abrupt transition around a critical connectivity k = 2. In 
this case the transition becomes sharper as the system size 
increases, suggesting a first-order transition type of behav- 
ior. The inset shows the total number of redundant bonds 
divided by the number of nodes in the network, N r = Nt/N. 
For a — 2.7, redundant bonds are not included before k = 2. 
In all simulations, we use /3 = 1.0 and take an average over 
1000 realizations of the growth process. 



all clusters have approximately the same size S. How- 
ever, without redundant bonds, the clusters are all tree- 
like with exactly 5—1 bonds, and t = 0. At this point, 
we have AH r = S a , and AH m = 2S 2 . Thus, for large S, 
the inclusion of merging bonds will lead to smaller energy 
variations, as long as a > 2. We then conclude that, in 
the large cluster limit, S> 1, both scenarios are stable 
for a > 2. 

The evolution of the system towards tree-like or fully 
connected clusters is determined at the beginning of the 
growth process. Considering that S = 3 represents the 
minimal size necessary for the inclusion of a redundant 
bond, we obtain AH r = 3 Q and AH m = 2 x 3 2 = 18. 
Thus, merging bonds become more probable when a > 
ln(18)/ln(3) = 2+ln(2)/ln(3) « 2.63, which corresponds 
to a threshold condition above which the system exhibits 
an abrupt transition. One should note that this is an 
approximate result, since we do not account for fluctua- 
tions in the cluster size distribution. However, as shown 
in Fig. [21 the results for a = 2.5 and 2.7 indeed con- 
firm the change in behavior from a sharp transition for 
the larger value of a to a slow continuous growth for the 
smaller value. Note also that the threshold value for a is 
not universal and could be readily changed by adding a 
multiplicative constant to any of the two terms constitut- 
ing the Hamiltonian of Eq. (JTJ> . In the inset of Fig. [21 we 
show the dependence of the fraction of redundant bonds 
N r = Nb/N on the average connectivity of the network 
k. As one can see, for a = 2.7, the inclusion of redundant 
bonds is delayed up to k ~ 2, confirming that the system 
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FIG. 3: Growth process when redundant bonds are favored. 
Here we show results for = 1.0 and a = 2.0. Since in this 
situation merging bonds are less likely to be included, the 
graph has to reach states where it splits in several fully con- 
nected sub-graphs, before a new merging bond is introduced. 
When the merging bond is included, a new and larger cluster 
is created. This explains the presence of discontinuous jumps 
in the size of the largest cluster. Assuming that the system 
consists of only fully connected clusters of the same size, we 
obtain the dotted line shown in the inset, S = fc + 1. This con- 
dition corresponds to the minimum bound for the simulation 
results, that approximately follows this theoretical prediction. 
Since the largest cluster S is finite for any finite connectivity 
k, the system does not display a percolation transition. 

is in the tree-like regime. 

Let us examine in more detail the scenario for a small 
value of a = 2.0. In Fig. [3] we show that the fraction oc- 
cupied by the largest cluster P^ systematically increases 
with the average connectivity k, with a growth rate that 
decreases with system size N. The inset of Fig. [3] shows 
the same results, but for the size of the largest cluster 
S = NPca. One can see that S follows approximately a 
linear growth with the connectivity k. In this scenario, 
a merging bond is expected to be placed only when all 
clusters become saturated with redundant bonds. If we 
now use that all clusters have about the same size S, we 
obtain k = 5 — 1, which corresponds to the dotted line in 
the inset. The deviations of the numerical results from 
this prediction should be expected. In the growth model, 
the merging of two clusters can only double the value of 
S, so that the values of S at the plateaus observed in the 
curves are approximately powers of two. However, we 
see that the curves always approach the dotted line be- 
fore doubling S. This linear growth for with a slope 
that decays with the system size N, indicates that, in 
the thermodynamic limit, this system does not undergo 
a percolation transition, namely, Poo — for any finite 
value of k. 

Figure [4] shows results for a = 3.0. Here we are in 
the scenario where the clusters grow as loopless trees. 
In this case, the system undergoes a transition that be- 
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FIG. 4: Growth process when merging bonds are favored. For 
/3 = 1.0 and a = 3.0, the system does not include redundant 
bonds and all clusters remain tree-like until the critical point 
k c — 2 is reached. Supposing that the system comprises only 
clusters of the same size S, we have that S = 2/(2 — k) for the 
case of trees. This relation works as a minimum bound for the 
simulation results, as shown in the inset on the left. Thus we 
have the critical condition k c — 2, where the size of the largest 
cluster diverges to occupy the whole network. For k < k c 
the largest cluster remains finite and its occupation fraction 
Poo goes to zero as the system size grows, characterizing a 
typical first-order transition. The inset on the right shows the 
threshold connectivity kt to obtain a largest cluster greater 
than the square root of the system size, S > N 1/l2 (black 
circles), and greater than half the system size, S > N/2 (red 
squares). The red line follows k = 2 — 4/iV, the expected 
behavior for the connectivity where S = N/2. The black line 
is a fit of the form k = p x + p 2 x N~ P3 , with p 1 = 1.99 ± 0.03, 
p 2 = 3.14 ± 0.02, and p 3 = 0.53 ± 0.05. In the limit N -» oo 
both curves converge to k f» 2, that is in the thermodynamic 
limit we observe at k = 2 a discontinuous transition in the 
order parameter from a vanishing fraction, Poo ~ iV" 1 / 2 , to 
a finite fraction, Poo = 1/2, confirming the approach to a 
first-order transition. 



comes sharper as the the number of vertices N increases. 
Again, if we assume that the system is divided in trees 
of the same size S, we obtain S = 2/(2 — k), as indicated 
by the dotted line in the inset. As before, the size S in- 
creases in steps due the out-of-equilibrium nature of the 
growth process. Strikingly, the theoretical relation for 
the equilibrium state still provides a consistent predic- 
tion for the lower bound of the largest cluster size. Since 
S remains finite for any k < 2, and at the critical point 
k — 2 a tree that spans all the system is formed, it fol- 
lows that the order parameter P^ displays a first-order 
transition in the thermodynamic limit. 

As already mentioned, this growth process bears some 
differences with a thermal equilibrium state of graphs 
with the proposed Hamiltonian Eq. (TTJ) at low temper- 
atures. In fact, for k — ► 2 there is always an energy 
gain in breaking large trees in smaller highly connected 
graphs. One may then ask whether the sharp transition 
observed in the simulations is just a feature of the ir- 
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reversible growth process or could be reproduced in an 
equilibrium statistical framework. We now show that in 
fact in the limit of large a we can obtain an exact equi- 
librium solution that exhibits a first-order transition. 

If we impose that all clusters in the system are loopless 
trees, a — + 00, it is possible to enumerate all possible 
ways in which the network can be divided in a set of 
clusters of a given size. Let f2 represent the number of 
ways that a fully connected graph can be divided in trees 
with rii trees of size i = 1, 2, 3... We then have 

where N is the total number of vertices in the network, 
and T s is the number of trees that span a fully connected 
graph of size s, given by Cayley's formula, T s = s s ~ 2 [T^ |. 
Since all clusters are trees, we can relate the number 
of clusters N c with the system size N and the average 
connectivity k as N c — N(l — fc/2). Therefore, given 
a fixed value of k, the values of n s obey the following 
two constraints: ^n s — N c and ^ sn s = N, where the 
sum is over all possible cluster sizes s. In our generalized 
percolation model, we still need to impose a fixed energy 
value, ^2E s n s — E, where E s = s 2 is the energy of a 
tree with size s. Using Lagrange multipliers, 77, A, and 
(3 to deal with each of these constraints, we can find the 
cluster size distribution that maximizes f2, 

n s = e"^ s2 - As ^. (3) 
si 

The critical condition takes place when the distribu- 
tion Eq. © diverges. One can verify that, for (3 = 0, 
this happens when A = A c (/3 = 0) = 1. The critical 
connectivity can then be determined as 

*-('-$)-'('-&)■ » 



yielding k c = k(X c = 1,(3 — Q,i] — » 00) = 1. Note 
that, at the critical condition n s ~ s~ 5 / 2 , the fraction 
occupied by the largest cluster follows = k — k c , thus 
reproducing the known critical properties of the standard 
Erdos-Renyi model [la ]. 

For (3 > 0, the distribution always converges unless 
A — > —00. From Eq. J3J), we obtain that k c = fc(A — » 
—oo,f3,ri — + 00) = 2. For k < k c all clusters are finite 
trees, therefore occupying a vanishing fraction of the net- 
work. At k = k c = 2, a giant tree spans the entire net- 
work, characterizing a first-order transition. Of course, 
this simple approach to the problem is only possible due 
to the imposition of tree- like clusters. The general enu- 
meration of connected graphs with any number of redun- 
dant bonds is not a simple task [16j |. and the cluster size 
distribution in this generalized condition might be quite 
different. However, at least in the situation where redun- 
dant bonds are not present, explosive percolation can be 
duly explained within the framework of equilibrium sta- 
tistical mechanics. 

In summary, we have shown that two simple condi- 
tions, namely, the absence of loops and the imposition of 
clusters of similar sizes, are the only necessary ingredi- 
ents for a percolation process to display first-order tran- 
sition in the size of the infinite cluster as function of the 
average degree of the network. We argue that both con- 
ditions arc implicitly present in the explosive percolation 
model proposed in Ref. [lof . We emphasize that these 
conditions are essentially non-local, namely, the proba- 
bility of adding a particular bond depends on the global 
structure of the graph. Moreover, our model provides 
a simple connection between explosive percolation and 
equilibrium statistical physics, leading to a clear inter- 
pretation of the mechanisms behind this growth process. 
Finally, other possibilities for the energy function can also 
be investigated in different contexts, revealing a whole 
new family of percolation- like models. 
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